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I— i' We introduce a measure of decoherence for a class of density operators. For Gaussian 

■ density operators in dimension one it coincides with an index used by Morikawa (1990). 
I , Spatial decoherence rates are derived for three large classes of the Galilean covariant 

I quantum semigroups introduced by Holevo. We also characterize the relaxation to a 

^ ' Gaussian state for these dynamics and give a theorem for the convergence of the Wigner 

. function to the probability distribution of the classical analog of the process. 

(N ■ 

^ : 1 Introduction 

en ■ One important phenomenon in quantum optics is the suppression of wave behavior for a quan- 
^ 1 . tum particle interacting with an environment. This subdued wave behavior is usually referred 
^ ! to as decoherence and is strongly emphasized by many physicists [20] as being a major ingredi- 
I ent for the construction of a macroscopic world that is well-approximated by models of localized 
^ I objects following well-defined trajectories. Apart from the natural theoretical appeal of this 

■ topic, quantifying spatial decoherence has also attracted interest from experimental physicists 
^ . working in quantum optics [T71 [TT]. 

If pt is the reduced density operator of a particle with spatial degrees of freedom at time t 
interacting with a environment, then the rough intuition is that the particle is undergoing spatial 
decoherence if the off-diagonal position ket entries xi ^ X2 of pt{xi,X2) vanish at exponential 
rates. Thus the particle decohering through an environmental interaction is in some sense 
becoming more diagonal in the x-basis. In the present paper, we study certain categories of 
dynamics for decoherence by introducing a coherence index of the form: 

CAP) {ik^Un-ix^^p? 

DAp) I 1 



SAp) = = ^^"^^ ' TT' (1-1) 
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E-=iTr[{A,-Tr[A,p],pP] 
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for a density operator p, where Xj j = 1, - ■ - d are the position operators for a particle traveling 
with d spatial degrees of freedom. The numerator CAp) is a coherence length-like quantity 
while the denominator DAp) is a standard deviation-like quantity. We study the above index 
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for a density operator Tt{p) in the limit t ^ oo, where Ft is a dynamical semigroup of trace 
preserving maps formally satisfying the equation: 

-^E4?[^-[^.'r(p)]]+ /"dMx,k)(w;v.(p)vr.,k-r,(p)). (1.2) 

In the above, K is the vector of momentum operators, Xj for j = 1, . . .d are the position 
operators, VF(x,k) = e'^^^'^^^^ is the Weyl operator corresponding to translation in phase space by 
{q,p), /X is a symmetric measure about the origin on M'^ x Mf^ satisfying J dp.{^, k)(|xp + |kp) < 
cxD, and y4^'^, (^A^'^Y = A^'^^A''^^ are the d x d block matrices of a positive semidefinite real 
valued matrix A: 



A 



The dynamics Fj describes a free particle (no forcefield potential) in a random environment 
giving the particle a Levy process of phase space kicks through conjugation by the Weyl op- 
erators. The quadratic terms in [Xj, ■] and [Kj, ■] correspond to a continuous limit of frequent 
small kicks. Later in this introduction, this model and related models will be discussed further. 
Define the 2o!x2(i matrix: B = J dfi{^,k) (k)®(k)- Let S^'^, S^'^, S*^'He the blocks of 

B: gl.k^. The analysis of the asymptotics of S^{Tt{p)) splits into three main categories. 

Let be some positive measure on M'^. 

1. Only jumps in momentum: A^'^ = A^'^ = A^'^ = and /u(x, k) = 5(x)zv(k), where A^'^ 
is assumed to be positive definite or ly is assumed to have a density. 

5,(r.(p))^r^V3 ^^(f + (1.3) 

2. Only jumps in position: A^'^ = A^''' = A^'^ = and //(x, k) = zy(x)5(k), where A''''' is 
assumed to be positive definite or u is assumed to have a density. 

%(r,(p)) ~ t-^2^Tr[(A^''=)-i]^^J e n i m^!^ ^^-^^ 

(J t^k|p(k,k)|2)2 

3. Active presence of both jumps in momentum and position where A is assumed to be 
positive definite or p is assumed to have a density. 

%(r.(p)) ^ r^VTi^^Mp.^ (1.5) 

For two functions at, Pt, by ~ (3t we mean that limj^oo % — expected that 

the asymptotics will have an error on the order of 0{t~^) for cases ([1]) and ([3]) and 0{t~^) 
for case ([2]) due to the application of variations of Laplace's method in the approximations. 
Similar statements can made for decoherence rates in momentum. 
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Theorem fl4.4p characterizes the relaxation for the dynamics Tt in the cases (1) and (3) 
above. It is essentially a central limit theorem where the Poisson noise can be approximated by 
a Gaussian noise and all information from the initial state is lost. For dimension one, case (1) 
and certain generalizations are discussed in [22] and [2^. In the current article, we characterize 
the relaxation as a Hilbert-Schmidt norm convergence of the process Tt{p) for an initial p to a 
process pt defined as 

Pi(xi,X2) = (^^^(^*))' ^-i(xi-X2|at(xi-X2))-i(xi-X2|bi(xi+X2))-i(xi+X2|ct(xi+X2))^ ^^ q^ 

(27r) 2 

where = t{A^'^ + B^'^), ht = |/, and q = ^{A^'^ + Since Tt{p) and pt tend to zero 

in II ■ II 2, in (14.41) we characterize the convergence for t — ^ oo as 

l|rt(p) - pth 



iPth 



0. (1.7) 



As discussed in [22] and [27], non-Gaussian limiting dynamics corresponding to more general 
stable laws can emerge if the Levy measure i/ in case (1) does not have finite second moments. 

Also for cases (1) and (3), Theorem (14. 5 p characterizes the convergence of the Wigner 
distributions Wrt(p)(x,v) to the probability distributions of the classical counterpart of the 
dynamics. Let {xt,kt) G M^" be the stochastic process in phase space with characteristic 
functions 

where 

/(q,p) = _l(p|A^.^p) _ (plA^-'^q) - l(p|A'=''=p) + J dp{x,k){e'P-''+''i-^ - 1). 

The stochastic process {xt + /q dsks, kt) describes a classical particle initially at the origin with 
zero momentum undergoing free evolution interrupted by phase space jumps determined by 
the Levy process (xj, kf). Let Pt(x, v) be the probability distribution at time t for the process 
{xt + Jq dskg, kf), then (14.51) states that 

(/^xc/v|Wr,(,)(x,v)-p,(x,v)p)^ 



( / (ix(iv|pt(x, v)|2) 



The index 5'^(p) is intended as a measure of the spatial coherence of the density operator 
p. In terms of the integral kernel of p in the X basis, S^{p) can be written: 



( / (ixi(ix2|xi - X2np(xi,X2)p 
S^{p) = ^ '—^ , for Tr[A:p] = m. (L9) 

^ / (ixi(ix2|xi + X2 - 2mP|p(xi,X2) 



The numerator measures a width of the integral kernel p(xi,X2) in the off-diagonal direction 
while the denominator measures the width in the diagonal direction. Since the diagonal of the 
kernel p(x, x) has the interpretation as the probability density of finding the particle at the 
point x, the computation of the denominator is centered by the mean value Tr[Xp]. S-g{p) 
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takes values in the interval (0, 1] and is one in the case that p is a pure state. For the case 
when p has a Gaussian form in one dimension 



^ - A{xi-X2)'^ -iB{xl~xl)-C{xi+X2Y -iD{xi-X2)- E{xi+X2)-F 



where all constants A, - ■ ■ ,F are real, A > C > 0, and F = then 5'^(p) = y ^- This agrees 
with the index used in [23] for decoherence models where the Gaussian form is preserved. A 
similar quantity coherence for momentum S^{p) can also be defined. 

We are interested in the power laws that arise in the study of S^{p) when p = pt evolves 
according to some Markovian dynamics. In the case where pt evolves as a free particle without 
noise (i.e. equation (11.21) with zero noise terms), then Sj^{pt) will approach a non-zero constant. 
We see from the difference between (II. 3p and (11.41) that the power law for Sji{pt) will depend 
on the nature of the noise. 

There are bounds involving the numerator of S^{p) and the denominator of S^{p) and vice 
versa, since using the identities 

p= -i \Xj ,{Kj, p}] + {Kj,i [Xj , p] } and p = 2 [Kj , {X^ , p}] - {Xj,i [Kj , p] } , then 



CAp)Dg{p) > \ and Q(p)Z^^(p) > \. (1.10) 

The inequalities (11.101) reduce to the uncertainty principle in the case where p is a pure state, 
and they represent a fundamental lower bound on the coherence distance in position and the 
statistical spread in momentum in terms of the spread in momentum and position respectively. 
The coherence length may only vanish as the spread in momentum grows without bound. For 
more realistic dynamics [25] including the effects of friction, the wavelength will tend to settle 
down to some finite value; the suppression of wave effects will depend on the wavelength being 
small compared to relevant length scales. 

A more general understanding of the interplay between coherence length and momentum 
distribution can be found by looking at expression for C^(p) = ||p||2 ^|| [-^5 p]||2- Since 

Vk(e^''^pe-^'^^)|k=o = ^[X,p], 

— * 

z[X,p] gives the rate of change of the density operator p given a small momentum shift in 
momentum k. C^{p) is thus susceptible to being small in situations where the momentum 
distribution of p is widely spread out and smooth. 

Unfortunately, the measure of the momentum spread appearing in (ll.lOp is not simply the 
standard deviation: (Tr[p(K — Tr[ii'p])^]) 2 or a quantity naturally bounded by it. We choose 
the quantity -D-^(p) for our coherence index S^{p), since it bounds Dj^{p) and satisfies (11.101) . 
We could have also considered alternative indices 

II[^,Vp]||2 , Il[^,v/p]ll2 

Y cLllQ Y 1 

2(Tr[p(X - 1i[Xp\Y]Y 2 (Tr [p(i? - Tr [i?p]) 2] )2 

which take on values in [0, 1] and have the same numerator-denominator inequalities as in (ll.lOp 
hold. The advantage of these quantities is that the denominator is actually a variance formula. 
However, the numerator appears difficult to work with. 
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The convergence fll.7p gives a detailed form of the hmiting dynamics. We can for instance 
write down the form of the density operator pt in the momentum basis: 

Pi(kl,k2) = (^^^(^f))' ^-^(kl^k2|a^(ki-k2))-A(ki-k2|b;(kl+k2))-^(ki+k2|c^(kl+k2))^ 

(27r) 2 

where = + B^'^), b[ = — |/, and c[ = jj.{A^'^ + B^'^)~^. Since the eigenvalues 

of c'f. are vanishing as t ^, the average kinetic energy of the particle is growing linearly with 
t. This corresponds to a high reservoir temperature assumption appearing in derivations of 
various subclasses of dynamics satisfying (11.21) . The off-diagonal width of pt(ki,k2) is on the 
order of t~2. With the naive identification ki ^ k2, (II. lip is like a Gibbs state A^^^^e^^'l^'^, 
where = (, where we have assumed that the environment is rotationally symmetric so 
that A^'^ + B^'^ = cl for some constant c). Thus the "temperature" increases linearly in time. 
As a consequence of the Gaussian form of pt 

so the particle decoheres in momentum at the same law that it decoheres in position. 

(II. 8p gives another way to define the classicality of the system as t ^ oo other than 
the vanishing of the off diagonal terms for the density matrix prt(p) in the position and the 
momentum basis. By taking the Fourier transform of pi(xi,X2) in the xi — X2 direction we 
can obtain the relaxation form for the Wigner function >Vrt(p)(x, v) and the classical process 
Pt(x, v). 

Special cases of the dynamics (11.20 have been derived in the study of decoherence by various 
authors. In [TFj , the authors discuss the reduced dynamics Tf for a spinless particle interacting 
with a gas under the assumptions that the reservoir of gas particles is translation invariant, 
interaction particles from the reservoir are in an ensemble of momentum states (commuting 
with the momentum operator), the reservoir is not effected by collisions with the particle, 
collisions are instantaneous, and an additional length scale assumption about the collisions the 
particle receives. In the three dimensional case, the derived Schrodinger dynamics take the 
form: 

|r,(p) = ^[\K\^Up)] - I [X„r,(p)]]. (1.12) 

The first term on the right is merely the free dynamics generator, but the second term on 
the right represents the stochasticity introduced by the reservoir. Equation (11.120 describes a 
free particle interrupted by Wiener motion of jumps in momentum. Notice that the generator 
has the Lindblad form with irreversible part: L{p) = XpX — \X'^p — \pX'^. Looking at 
operator elements in the x-basis, L{p){xi,X2) = — ^(a^i — X2)^p(xi, X2), so the stochastic term 
indeed seems to generate an exponential vanishing of off diagonal entries. Intuitively this effect, 
however, is somewhat mitigated by spreading out from the free dynamical term. An analysis of 
the decoherence of this model in dimension one is studied in [19] and also in [23], where some 
additional terms in the Lindblad form corresponding to a harmonic oscillator potential and a 
friction term are also considered. In the analysis of [13 [23], it assumed that the initial density 
operator p has a Gaussian form: 

^ -A{xi-X2)'^-iB{x'(~x'^)-C{xi+X2)^-iD{xi-X2)-E{xi+X2)-F 



5 



where all constants A, - ■ ■ , F are real, A > C > 0, and F = j^. is the standard deviation 

of the Gaussian state in the position variable. The quantity is interpreted as the coherence 
length of the state. In quantum optics, the coherence length is the approximate length at which 
different parts of the wave packet interfere. The authors in fL9\ use the fact that the dynamics 
Tt maps Gaussian density operators to Gaussian density operators and derives differential 
equations for the coefficients At, ■ ■ ■ ,Ft. The relevant quantity for the study of decoherence 

is the asymptotics of the ratio (which is equal to the coherence length divided by the 

standard deviation at time t). For the model fll.l2p . the asymptotics are ~ ct~^ for some 
constant c. 

In [9] there is derivation closely related to that in [19], but without the short length scale 
assumption. The derived dynamics satisfy a differential equation which can be written 

jTtip) = z[\K\', F,(p)] + J n(k)rfk(e^'^^F,(p)e-'^^ - F,(p)), (1.13) 

where n{k) is a positive density. These dynamics describe a free particle with a Poisson field 
of jumps in momentum, where jumps by k in occur with rate ?T,(k)(ik. An equation of this 
form was originally introduced in [10] as a fundamental alternative to the Schrodinger equa- 
tion rather than an effective reduced dynamics for a particle interacting with an environment. 
The dynamics have also have been used to make quantified comparisons with the results of 
experiments [T6| IT]. For a general discussion of decoherence with an emphasis on these models 
see [2D|. 

The dynamics described by (11.121) and fll.l3p both share the property that they correspond 
to an environment that is homogenous. In fact they both satisfy the covariance relation 

x+tk,k), 

(1.14) 

for all Weyl operators W(^x,k)- This follows because conjugation by W^^k, which corresponds to 
shift in phase space, commutes with the noise part of the generators. Moreover, if Ft is the free 
evolution generated by i[K^,-], then Ft{Wl^^,^^pW(^x^k)) = W*^^^i^^^-^Ft{p)W(a:+tk,k)- Hence, even 
after time evolution, conjugation by Weyl operators commutes with the noise. A Schrodinger 
dynamics F^ satisfying (11.141) is said to be the Galilean covariant. 

Intuitively, a Galilean covariant semigroup corresponds to a free particle traveling in a 
random environment that is invariant with respect to translations in phase space. In other 
words, the probability of the particle undergoing a sudden shift (A(x), A(k)) in its position 
and momentum is independent of its location in phase space. In [12j, there is a complete 
characterization of these processes in terms of their Lindblad form with an additional assump- 
tion that the dynamics satisfies rotational covariance Tt{U*pUa) = U*Tt{p)Ua, where cr g 5*03 
and ([/o-/)(x) = /((Tx). Although Holevo worked in the Heisenberg representation, in the 
Schrodinger representation the dynamics formally satisfy: 

jTtip) = ^[\K\^Tt{p)] - \ X^(c^'^[X„ [X,,F,(p)]] +c^'nx„ [i^„F,(p)]] 

+ c^^^[K,,[X,Jtm+ c^''[K,,[K,Jtm) + j dp{^,'k)[W^,^Tt{p)W:,y-Tt{p)l (1.15) 
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where the matrix (^k,x '^k,k^ is real-valued and semi-positive definite, the measure is on 

X M'^ — {0} has the rotational invariance //(x, k) = /i(crx, ak), and the measure has the Levy 
condition: 

/ ^^^^''^ ^ + |Ip + |kp ^""- ^^-^^^ 

For the integration in fll.lSp to make sense, the integration is taken over spheres centered at 
the origin first and then in the radial direction (to get a quadratic weight from the integrand 
near zero). Looking at the application of Proposition (2) in fT2] to the proof of the Theorem 
on page 1819 of [12j, we can see that the analogous results hold when the rotational invariance 
is removed. In the case when the rotational symmetry is replaced by just the origin symmetry 
(i^(x, k) = (i/i(— X, — k), the corresponding dynamics can be written in the form (11. 2p . 

The case in which there is only a Poisson term and fj,{x,k) = 6{x)i'{k) corresponds to the 
form derived in |9], and the case in which there is no Poisson term and only the <f'^ quadratic 
term is non-zero is the model derived in [19]. Larger classes of dynamics which notably do 
not include the energy blow up implied by the form (II. lip are derived in [25] from scattering 
arguments. For a survey of various dynamical semigroups corresponding to spatially invariant 
environments, see p6J. 

One interesting aspect of Galilean covariant dynamics is their constructibility using clas- 
sical stationary stochastic processes with independent increments. The dynamics Fj can be 
constructed as 

r*(p) = (1.17) 

where Vt is the unitary group is generated by |Kp (free dynamics) and {xt, kt) is a stationary 
stochastic Levy process with characteristic functions. 

where 

/(q, p) = _lc^.^|p|2 _ c^''=q . p - ic'=''=|p|2 + j dti{x, k)(e^P-''+^'i-'^ - 1). 

Note that we have stated the result (11.171) for the dynamics Tt, but the construction in [I2] was 
made for the adjoint dynamics FJ' = (Heisenberg representation). The dynamics Tt are thus 
a statistical average over certain unitary trajectories constructed using the Weyl operators and 
the free unitary dynamics Uf. 

The closed factorized form of the characteristic function as found in (I1.19P for the quantum 
characteristic function of the covariant dynamics is a consequence of the constructibility of 
the dynamics Ft using only conjugation by Weyl operators and the Vj's. It is shown in [13] 
that this implies that Weyl operators evolved under the adjoint dynamics and can be explicitly 
computed as: 

$t(^(q,p)) = e^o 'i-'(q+(t-)p.p)i^(q^,p^p). (1.18) 

A discussion of the dilation of the full collection of processes described by Equation (ll.lSp can 
be found in For a larger discussion of dilation of quantum semigroups using classical noise 
see [T3] . 
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In our analysis, we make use of the fact that the map of a trace class operator p to its 
quantum characteristic function ¥?p(q, p) = Tr[VF(q_p)p] extends to an isometry of Hilbert- 
Schmidt operators to L^(M'^ x M'^, j^^^dc]dY)). The dynamics formally satisfying (11.21) has a 
closed expression for the time-evolved quantum characteristic functions '^Vtip)'- 

¥'r,(p)(q,p) = L2\V p j\ \ p j / ^.^^ ^ ' 'Vp(q + ip,p), (1-19) 

with -0^(q, p) = y" (i/i(x, k) (cos(q • k + p ■ x) — l) . (1.20) 

It is shown that the expression ^^(q, p) can be effectively replaced for the sake of computing 
the asympototics of S-^{Tt{p)) with the quadratic form from the second-order Taylor expansion 
of '0/i(Q)P) the origin (q, p) = (0,0). This approximation has the flavor of a central limit 
theorem where the Poissonian noise can be replaced by a Gaussian noise and can be seen 
in (14.41) . Laplace's method makes it possible to find the asymptotics of quantities needed to 
calculate Sji{Vt{,p))- 

Note that the equation (11.191) is identical to the characteristic function for the stochastic 
process {xt + /q dskg-, h) except that v?p(q, p) should be replaced by the characteristic function 
for (xo, ko) (in the general case, although we have assumed that the process began as 6 dis- 
tributed at the phase space origin). The non-classicality arises through the initial characteristic 
function v'p(q, p) which will not have the form of a characteristic function for a classical prob- 
ability distribution over phase space. However, as t ^ cxo, the multiplication factor in (I1.19p 
vanishes away from origin where V'p(q, p) takes the generic value of one. This is the basic 
reasoning behind the convergence (11.81) . 

This article is organized as follows: Section [2] gives a general discussion of coherence indices 
of the type 11.11 Section [3] gives a brief background on the meaning behind the formal Lindblad 
equations with unbounded generators as found in the work of Holevo [121 [IS El IIS]) cind 
Section H] contains the main results of this article. 



2 State Coherence Indices 

Let 7i be a complex Hilbert space and Aj, j = 1, . . .dhe a. family of self-adjoint operators with 
essential domains Dj, and let p be a density operator such that Ajp is trace class {Ajp G Ti(7i)). 
Define C(^Aj)ip) and D(^^.)(p) through the following formulas 

C^(A.)(P) = ^(E||[^.,P]||0' (2.1) 

and 

D^aMp) = ^(EllH-Tr[A,p],p}||^)'. (2.2) 

The operator Ajp is defined through the bounded bilinear form B{g,f) = {Ajg\pf). Notice 
that p maps arbitrary vectors / to the domain of Aj. 

D{Aj){p) is intended as a sort of standard deviation for the operators Aj in the state p, while 
C(Aj){p) gives some sort of measure of how close the family of observables Aj are to commuting 
with the state p. 
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Definition 2.1. Let Aj be self-adjoint operator with dense domains Dj and p G Ti(7i) be a 
state such that Ajp G T(7i) for each j. If D(^a ){p) 7^ 0, then the index S(^a ){p) of the family 
{Aj) with respect to the state p is defined as 

If tlie observables Ai liave some form of units (e.g. lengtli, energy), tlien the index yields 
a dimensionless parameter related to the commutativity of the observables A^ with respect to 
the state p. The following proposition gives a few basic facts about the index S(^Ai){p)- 

Proposition 2.2. Let Ti be a Hilbert space, p be a density operator, and (Aj), j = 1, . . . d be 
a family of self-adjoint operator with domains Dj such that Ajp is trace class. 

1. 5(^^,)(p)G [0,1]. 

2. S(^Aj){p) = iff P commutes with every Aj. 

3. If p = \ f){f\ is pure and f is not eigenstate of Aj for all j = 1, ■ ■ ■ d, then S(^Aj){p) = 1- 
Example 2.3. For a 1-dimensional Gaussian state: 

_ "^a/C^ A{xi-X2f-iB{xl-xl)-C{xi+X2)^-iD{xi-X2)-E{xi+X2)-F 
P ~ /— ^ ) 



Cxip) = Dx{p) = o.iT'd the quantity Sx{p) = agrees with the index used 

in IM\M- 



Although Cx{p) and Dx{p) differes from the quantities and interpreted as the 

coherence length and the standard deviation in [23] by a factor of 2~5 for Gaussian density 
operators, Cx{p) are Dx{p) are not amenable to an interpretation of this sort for a general 
state p. The squaring of an expression involving p as found in the formula for S^{p) can give 
a skewed weight for the probability weights of events. The following example gives an extreme 
situation where this becomes apparent. 

Example 2.4. LetH = L\R), = (/^^ Idx^hA^ where A.^ = EZi ^ EZi 

Define the density operators p„ = Y.m^n,m\4>m){(/>m\, with A„,i = ^, Xn,m = ^ for 2 < m < 
V? — v? + 1 and \n,m = otherwise. We can calculate the numerator of Sx{pn) using the 
following. 



\[X.Pn]\\l = ^m\X'^,) - (0l|X0i)2) + 0(4). 



Moreover we can approximate the denominator of Sx{pn) as. 



\\{X.Pn}\\l = ^{h\{X - Ifh) + - + 0(- 



Hence, 

((0l|X20i)-(0i|X0i)^ 



Sx{Pn 



((0i|(X-l)20i) + (0i|(X-l)0i)2). 



The above expression depends only on the first state (pi even though this state has a weight of 
only which is a diminishing fraction of the total weight. 
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In general, just as for classical diffusion processes, only states of very specific forms can 
occur after they have been acted upon by an irreversible environment p ^t{p)- The states 
that are likely to occur depend on the nature of the environment. Our analysis, in Section H] 
essentially relies on the fact that when stochastic shifts in momentum are present, then after 
sufficient time Tt{p) becomes essentially Gaussian-which is the statement of Theorem ( 14.4p . 
Thus for those dynamics, S^{p) is asymptotically expected to serve well as a coherence index. 
In the case where there are only stochastic shifts in position, the state Tt{p) becomes in some 
sense only partially Gaussian since the asymptotic characteristic function fll.lQp will only be 
forced to be quadratic in the exponent with respect to the p variables. This seems apparent in 
the asymptotics (11.40 . since the constant 

{jdk\pik,k)\'\^-E[Kp]\'y^ 

(/rfklp(k,k)p)l 

has the strange squaring of p(k, k). With a more accurate formula for the coherence length 
divided by the standard deviation in position, we expect that this constant would be replaced 
by a variance formula for the probability density p(k, k). In general, for other types dynamics 
we expect Sji{pt) to capture the correct power law even if the constants are not completely 
natural. 

In the proposition below, we give useful expressions for the quantities need to compute 
Sxip) and S^{p). 

Proposition 2.5. Let p he a state such that Jp G Ti(L^(]R'^)) for any 
J G {Xi, ■ ■ ■ , Xd, Ki, - ■ ■ , Kd}, and define 

Vp = (Vp(/Jp)(0, 0) and Vq = (Vq(/?p)(0, 0). 

Then for the vector of position observables X, 

l|[^,P]||2= ( J c/qdp|q|Vp(q,P)r)",«n(^ ||{X,p}||2= ( J (iqdp| (Vp - Vp)^p(q, p) ^ 

For the momentum variable K, 

\\[K.p\h = {j rfqf/p|pr|y.p(q,p)r)Vnd \\{K , p}h = { j rfqrfp| (Vq - Vq)^,(q, p) 1^) 



2\2 



Proof. By ( ]A.3[) . quantum characteristic functions define an isometry from Hilbert-Schmidt 
class operators to functions in L?'{W^ x W^, j^^^dpda^ (Lebesgue measure on phase space mul- 
tiplied by a factor j^^)- Hence we have that 



Tr[-[^,p]^] = j dqdpJ^lVilx^pM^P) 



2 



By definition ip^^j^ ^](q, p) = Tr[e*('''^+P'^)i[X, p]]. However, we can write i[X, p] = Vd\d=oWipWg: 
where convergence for the limits jiW^^^pWhe, — p) — ^ i[Xi,p] takes place in the trace norm by 
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Lemma (lA.ll) . Since the convergence is in the trace norm it follows that we can commute the 
limit with the trace in the following computation: 

Tr[e^('i-^+P-^)^[X,p]] = Vda=oTr[e^('i-^+P-^)l^(to)PW^(a,o)] 

= V(,,o)|a=oTr[iy(,^,o)e^^''-'^+P-^^W-(s,o)P] = V,|a=oTr[e^«'i-(^-'^)+P-^V] 

Hence we can conclude that 

/I 
dqrfp|q^,(q,p)p)'. 

The quatities ||{X,p}||2, ||[i?,p]||2, and [[{K, p}||2 have similar arguments. 

□ 

3 Covariant Quantum Dynamical Semigroups 

For the purposes of the decoherence analysis in the next section, we work with the character- 
istic functions fll.lQp . using fl2.5p . In this section, we discuss the meaning behind the formal 
Markovian master equations with unbounded generators discussed in the introduction. For a 
more in depth view of this topic see fl2[ [T3t [HI [15] and further references. We finish up by 
making a few comments on the action of the dynamics from the perspective of characteristic 
functions. 

Given an Hilbert space 7i, a dynamics can be seen as a collection of completely positive 
maps (cpm's) in the Schrodinger picture acting on trace class operators : Ti(H) Ti(H), 
or in the Heisenberg picture acting on bounded operators $t : B{T-C) B{7i). The dynamics 
Tt and $t are related through the trace formula: 

Tr[F,(p)G]=Tr[p<|.,(G)]. (3.1) 

Since T(7^)* = B(H), the maps Ft are pre-adjoint to Although physicists working in quan- 
tum optics tend to work in the Schrodinger picture, those working on existence and uniqueness 
of Lindblad type equations tend to use the adjoint dynamics. Through Equation (13. ip either 
dynamics can be constructed using the other. Also, the dynamics F^ is trace preserving iff 
the adjoint dynamcis $j is unital (i.e. $t(/) = / for all t). The maps $j are said to form a 
dynamical semigroup if $t$s = ^t+s, and Tr[p$t(G)] is continuous (i.e. weak*-continuous). 

In p[4j, Holevo studies a dynamics $t operating on B{L'^(R^)) in the Heisenberg represen- 
tation and formally satisfying: 

j 

-c^'^[K„[X„$t(G)]]-c'^'^[K,[K,$t(G)]]) + J dfii^,k)[W^,^MG)Wl^- MG)], (3.2) 

^k,x ^k,k j is a positive matrix with real valued entries and is a measure on x R-^ 

satisfying the Levy condition / dn{x, k) i_^xp+jkp ^ ^ rotational invariance /u(x, k) = 

fj>{cr:s., crk) for a G 5*03. 
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Since the Lindblad Equation (13.21) has an unbounded generator, the classic result [21] guar- 
anteeing the existence and uniqueness of a norm continuous adjoint semigroup $f of completely 
positive maps satisfying = / does not apply. Just as in the case of generators of uni- 

tary groups, unbounded generators of Markovian semi-groups require extra care to define and 
pose new technical difficulties. One approach for dealing with these technical issues is the 
introduction of a form generator. 

Definition 3.1. Let D dTi he dense. A form generator is a linear map L : Dx B{li.) x D ^ C 
such that for f,gdD and G G B(H), 

1. 

C{g;G; f) = Cif;G*;g) 

2. 

J2 ^fi; GIG,- /,) > when ^ G,f, = 

1,3 3 

3. For any fixed g, f , C{g; G; f) is continuous in G with respect to the strong topology over 
any bounded subset of B{7i). 

A form generator C is said to be unital if C{g; I;g) = for all f,g & D. The definition 
for the form generators is inspired by the form of a bounded Lindblad generator. In [i2\, it 
is shown that for any form generator £ there exist operators Lj, j G N and B with domains 
including D such that 

C{g; G- f) = J2{Lj9\GL,f) - ^iBglGf) - ^{glGBf) 
j 

Given a form generator C, we can then ask if there is a process $t satisfying $o(G') = G 
and 

j^{g\MG)f) = C{g;G;f), (3.3) 

where g,f E D and G G B(H) and some regularity properties are assumed for Tt{G). An 
important criterion used for the construction of solutions to this equation is that i? is a maximal 
accretive operator. By an analogous result to Stone's Theorem [24j, maximal accretive operators 
are the generators of strongly continuous semigroups of contractive maps [18]. In [12] it is 
shown that for any unital form generator C admitting a Lindblad form where the operator B is 
maximal accretive then there exists a unique minimal dynamical semigroup $ to the equation 

f^{g\<!>,{G)f)=C{g;MGy,f) (3.4) 

where ^o{G) = G. A solution to the above equation is said to be minimal, if for any other 
solution 

^'tiG) > <^t{G) when G > 0. 

Surprisingly, the conservativity of the minimal solution ($*(/) = /) is not guaranteed if the 
form generator in unital. A general set of necessary and sufficient conditions for guaranteeing 
conservativity is unknown, and in the literature stringent conditions are assumed in order to 
prove the conservativity for a specific class of form generators [6l [12] . 
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For f,gED = n^^pDom{p - K + q- X), the form generator C{g] G; /) of the adjoint djTiamics 
corresponding to the formal Equation fl3.2p has the form: 

Cig; G- f) = Trig- G; f) + T2ig; G; f) + T.ig; G; /), (3.5) 

where 

T,{g; G; /) = -z{K'g\Gf) + i{g\GK'f), 

3 

T^ig; G;f) = J2 {c'^'^X^glGX,/) + c^''{X,g\GK,f) + c''^{K,g\GX,f) 

+ c'''''{K,g\GKjf) - ^((c"'"X| + c^'^iT^X,- + c"'%K, + c'''''K])g\Gf) 

- ^{glGic'^'^X] + c'^'^^K^Xj + c"'%ir,- + c''''' K]) f)) , 

Ug-.G; f) = J dfi{^,k){{Wl^g\GWl^f) ~ {g\Gf)), 

and the integral is taken over surfaces of equal radius to make the integration well defined. 
In [T^, it is shown that for an G G B{T-C) there is a unique conservative dynamical semigroup 
with ^o{G) = G and satisfying the equation 

j^{g\UG)f) = Cig;MG)f), 
and the dynamics have the covariance relations: 

**(W^(;,p)G'K/(q,p)) = W^(;+*p,p)'^'*(G')W^(q+*p,p), and UKGR.) = KMG)K. 

Conversely, it is shown that any conservative dynamical semigroup satisfying the covariance 
relations above is the unique solution to an equation of the form (13. 4p . 

Since the dynamics $i acting on any Weyl operator VFq^p is explicitly computable (11.181) . 
this imphes that the quantum characteristic functions of the predual process Tt are explicitly 
computable, since 

^r^p) = Tr[Wq,pr,(p)] = Tt[MW^,^)p] 

= e/o'(9+(*-^)p.p)'^^Tr[Wq+tp,pp] = e^o'('?+(*-«W)'^-<^^(q + tp,p). 

However, for the free dynamics Ft generated by 'fiFt{p){fl,p) = 'fpi^. + 'tPiP), hence 

in the formula above we have factorization of the quantum characteristic function with a noise 
part and a deterministic part. The stochastic factor e-^o K<i+{t-s)p,p)ds jg consequence of an 
analogous construction to (I1.17P for the adjoint dynamics $t and basic computations with 
Weyl operators. 

It is useful to think about how the dynamics act in terms of their quantum characteristic 
functions. We can define the action of the dynamics Ft and J-'t acting on characteristic functions 
through the formula: 

Av^p = ^Ftip), and J^t'^p = ipFtp- 
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Notice that Ft forms a semigroup of contractive maps on L^(M'^ x M.'^). This can be seen through 
the formula fll.lQp . but follows from more general considerations. The quantum characteristic 
functions define an isometry from T2(-L^(M'^)) to L^(M'^ x W^, (27r)~'^(ix(ik). The maps are 
completely positive and extend to unital maps (since they satisfy the same equation as the 
adjoint maps with i replaced by —i in the kinetic part of the Lindblad generator), so 

Taking the trace of both sides and using the isometry 

In many cases, we will find it convenient to write 

</'r*(p)(q, P) = Ftrl^piq, p), (3.6) 
where is the multiplication operator of the form 

r;^ = g/o'^^-KO p")l^("7'')>+'^-('i-^P'P)] . (3.7) 
4 Decoherence Rates for Covariant Dynamics 

In the section, we will compute decoherence rates for cases of covariant dynamics: where there 
is only stochastic shifts in momentum, only stochastic shift is position, and an active presence 
of both stochastic shift in momentum and position. The analysis is not an exaustive case 
analysis, since we always make an assumption such that either A^'^^^A^'^^A is positive (rather 
than just positive semidefinite) or that the measures p or u have densities. However, for the 
situation where p is assumed to have second moments, the main situations are covered. The 
proof of (14. 3p is essentially what is needed for (14. 4p . 

By the characteristic function isometry. Equation ( 13.71) . and the fact that Ft acts as an 
isometry on L^(M'^,M'^), we have that 

ll[^,r,(p)]||2 = -^||(q-tp)r;(^,)||2, and (4.1) 

[271)2 

\\{X - Tr[Xp], r,(p)}||2 = -^m^ + Vp - Vp^r.(p)(0, 0))r;(^,)||2. (4.2) 

(27r)2 

Moreover, by the origin symmetry of the noise, the noise does not change the expectation 
of the momentum and the position operators from the initial state. Hence with Ft = FtF^, 
E[Xp] = Vpipp{0,0) and E[Kp] = V,(/Pp(0,0), 

^p¥>F,riip) (0,0) = tVqV^r,p(0, 0) + V^y^rup) (0,0) = tVq(/^p(0, 0) + Vpcp.iO, 0) . (4.3) 
The last term from Equation (14.21) is bounded from above and below by, 

-^11 [tVq + Vp, r;]vpp||2 ± (-^11 (tVqv^p(o, 0) + v,¥.,(o, o))r;^ 

[2%) 2 ^[271)2 

+ -i-,||r,(tVq + Vp)v^,|h) (4.4) 
(27r) 2 ^ 

Where for some cases the later term will be seen to be of smaller order. 
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Proposition 4.1. Let p be a density operator such that J p G Ti(L^(M'^)) for all 

J E {Xi, ■ ■ ■ , Xd, Ki, - ■ ■ , Kd}. In the case when ^tAp) satisfies Equation M.19\} with p = 0, 



1. A'^'^ = A^''^ = A^'^ = 0, and A?^ positive, then 



Tr[y4^ 

2. A^'^ = A^'^ = A^'^ = 0, and A^'^ positive, then 

(/t/k|p(k,k)P)4 

3. and A is positive, then 



^ Tr[A^-j2 

Proof. For the numerator of S^(Tt{p)), we can use Equation (14. ip 

II [X, r,(p)] h = II (q - tp)rliip,) h, where 

[2nj 2 

r.'(^.) = pip/ dqdpe-/.'-<(-V')K-V)>-|^^(q,p)P. 

Computing the integral in the exponent and rearranging, 



I ds{(-V') \A {%■')) ds = i ((;) 1^ G)> + 1 (('if) \A 
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Now we will begin case analysis. 
Case 1: Since A^'^ is positive definite, there exists a unitary U and a diagonal D such that 
_ jj*DU gy changing variables t/q — >• q and f/(p — ^q) — > p, we can then write 

ll[^,r,(p)]||^ = ^ j dqrfp(i|qp + t2|pp)e-|('^l^'i)-¥(PPP). 

By Lemma (1A.2I) . ipp is uniformly continuous with (/?p(0,0) = 1. Hence if Aj are the entries of 
D we can apply Laplace's method to calculate the asymptotics of the above expression as 

[^(^ + ■ ■ ■ ^) + Jj^X^ + ■ ■ ■ ^)]ll^*('^p)ll2' Where ~ ^~''(Ai • ■ ■ A,)-^(^)'^. 

Hence \\[X Jt{p)]h t-"2i:T[A'^]h\\Tt{p)h. 

To get a hold of the denominator we study the expression -\ \\ [Vp — tVq, r'/]99p||2 from (14. 4p 

since the other terms have smaller order. The commutation is between derivatives and a 
multiplication operator and hence can be exphcitly computed: 



{2'kY J 3 
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If we rewrite t^q— yp = yQ— y(P ~ ^^)' then making the same change of variables as above, 
Laplace's method gives 



The other terms from -^\\{tV^ipp{0,0) + Vp¥?p(0, 0))r/((^p) Ha and ^rH/^/l^Vq + 

Vp)v2p(q, p) II25 can be at most of order t. For the second term we use that the derivatives of tpp 
are continuous and uniformly bounded by Lemma flA.2p in order to apply Laplace's method. 

Hence SAUp)) ~ r^v/S^^^M^^. 
Case 2: 

Let A'''^ = U*DU where D is diagonal with entries Xj. By changing variables t/p — > p, we 
obtain the expression 

^ I rfqrfp(|qr-2tqp + t2|pr)e-*<Pl^P>|y.,(q,[/*p)|l 

However, the third term dominates since |p| ~ In the limit, t ^ 00, the exponential factor 
places a weight on the surface p = 0. By Lemma ( ]A.2p . (fp is uniformly continuous and by 
Laplace's method we obtain the asymptotic expression 

m,Upml-t'^i . . )-^\\n{p)\\l, where \\n{p)\\l = ^ \ I dc\^p{^)\\ 

t Ai H Ad (27r)2 ^ (Ai---Ad)2 J 

Moreover, ipp{q) = Tr[e^'i-^p] = / dke''^-^p{k,k), so J dq\ipp{q)\^ = J dk|p(k,k)p. 

For the denominator we need to compute — ^-||(tVq+Vp— tVqy7p(0, 0)— VpV?p(0, 0))r^'((/7p)||2. 

(27r)^ 

The term tVq commutes with r^{ipp). The terms Vp and VpV?p(0,0) will be of lower order. 
Although when Vp acts on r^'((/?p)(0, 0)), it brings down a factor of tp, |p| ~ t^^. 



HKtVq -tVq(^p(0,0))r;(y,,)||2 - ^/ (iqrfpe-*<Pl^P>|Vq^p(q,p) - (Vq^)p(0,0)y.pP. 



f2 

Again the Gaussian weight is on the surface p = 0. In the limit t —>■ 00 this is asymptotic to 



tH-2 



Vq(/Pp(q, 0) - Vq</.p(0, 0)vPp(q, o)!^. 



(27r)f det(A'=''=)^ 

The integral on the right can be rewritten as j (ik|kp(k, k) — E[K p\p{k, k)p. 
Case 3: 

Since A is positive definite, just the first term alone yields exponential decay for phase space 
points (q, p) away from the origin. By changing variables tp — > p, 

dqdpe-*««'-^(ft)>^»(("*^')K°*l'))|^,(q,lp)p. 



In the limit t 00, the contributions from terms including jp become negligible for the 
asymptotics. The multiplication factor of j on the variable p in |(/?p(q, |p)P can only make the 
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function more amenable to Laplace methods since the function is effectively spreading out in 
the p variable. Neglecting these terms gives the same asymptotics as the first case. 

□ 

In the next theorem we consider the case where there is also noise is also a Poisson contri- 
bution to the noise. First we have the following lemma about classical characteristic functions. 
For an origin symmetric measure positive v on M'^ satisfying J (iz/(k) < -^g define the 
function 

v>,(i) = y dz/(k)(cos(k-i)-i). 

when z/(M'^) < oo, V^^(l) = ^vi\) - vi^'^). 

Lemma 4.2. Let v he a 'positive and possibly infinite measure on such that v is symmetric 
about the origin, and 

j dv{\i)\k\^ < oo. 

Then the first and second derivatives are bounded and continuous and an absolute maximum 
occurs at the origin. Moreover, if B is the matrix of second moments B = J (iz/(k)k ® k, then 
for any e there exits a 6 such that for all \l\ < S 



(1151) - -|1|^ < ^,(1) < -^-^(1|S1) + |1| 



2 " ' 2' ' - ^ ' ' - 2 " ' 2 

// in addition v has a density, then the absolute maximum of ip^ is obtained only at the 
origin and for any e there exists a 6 such that for all \l\ < S, 

+ '\l|51) < V'.(l) < -^^(1|M). 



2 

(l-e) 



Finally sup|i|=5 -^^(1|51) > sup|i|>^ ^^^(1) . 
Proof. We can rewrite the expression for ipiy as: 



■ 2 /'cos(k ■ 1) — 1 



M^) = J dM\K[-^^ — j- (4.5) 

The first and second derivatives in 1 of the family of functions /k(l), where 

cos(k ■ 1) - 1 

/k(l) = ^2 ' 

are continuous and uniformly bounded. By our assumption on u, the measure defined by 
(iz/(k)|kp has finite total mass. It follows that tpi, is bounded with bounded and continuous 
first and second derivatives. -0,^ is real, centrally symmetric, and the first derivatives of ifji, are 
zero at the origin. Thus (pi, obtains a absolute maximum at the origin. If u has a density ^ 
then the absolute maximum is unique since this is the only time point 1 at which all the phases 
in the integral ( 14. 5p are aligned. 

D'^'i/j^{l) can be expressed according to the formula: 



5)2^^(1) = - / dp{k)k(g)ke 
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For 1 = 0, this expression is equal to —B. D'^tpiylV) is continuous in the operator norm since its 
components are continuous and all norms are equivalent over finite dimensional spaces. 
In the direction 1, we can write the second order Taylor expansion: 



V'.(l) =^^^(0) + VV^.(0)1+ ds dr{\\D^t/j,{r\)\) = ds dr{\\D^Mrm 

Jo Jo Jo Jo 

Since —D'^i/j^iJ) is continuous with respect to the operator norm and positive semidefinite at 
zero, it follows for any e there exists a 6 such that 

(1 - e)B - eh < -D^iP,{\) < (1 + t)B + eh 

for all |1| < e. Applying this inequality to the formula flA.2l) . we have 

-^(1 + mB\) - ^|1P < ^.(1) < -^(1 - e){\\B\) + '-\\\\ 

In the case where v has a density, then the matrix B is positive definite since the integration 
of terms cannot have its support over some lower dimensional space. By continuity of 
D'^ipijiX), for any e we can pick a 5 such that 

-l(l + e)(l|Sl)<^.(l)<-i(l-6)(l|51). 

Furthermore, we can choose a 5 small enough such that '0i/(l) is concave down for all |1| < 5 
(and hence decreasing radially from the origin), and such that any local maximum that is not 
the origin is less than inf|i|<5 Hence 

sup = sup ipu{\) < sup --^(1 - e)(l|Sl). 

|1|>5 |1|=5 \\\=5 2 

□ 

For the Levy measure define 



B = j d/i(x, k) ® (Z) 



Theorem 4.3. Let p be a density operator such that J p G Ti(L^(R'^)) for all 
J G {Xi, ■ ■ • , Xd, Ki, - ■ ■ , Kii}. Let (prtip) satisfy Equation U.iy\) . where fi has origin symme- 
try /i(x, k) = X, — k) and the second moments 

J (iz^(x,k)(|x|2+ Ikp) < OO. 

1. With fi = 5(x)z/(k) and A'''^ = A^''^ = A^'^ = 0, then if v has a density or A^'^ is 
positive definite, then we have the same asymptotics as (1) from ^4-1^ with A^'^ replaced 
by A^'^ + S^'^. 

2. With fi = u(x)S(k) and A^'^ = A'^'^ = A^''^ = 0, then if u has a density or A'^''^ is 
positive definite, then we have the same asymptotics as (2) from with A^'^ replaced 
by A^'^ + S^'^. 
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3. If has a density or A is positive definite, then we have the same asymptotics as (3) 
from d^. j[ j with A replaced by A + B. 

Proof. The basic idea of this proof is that for long time periods we can effectively approximate 
the exponent of the expression ds'tp^{q^ — sp,p) as a quadratic through a Taylor expansion 

of iprnu at zero. Once we have shown this, then we can refer to our results from proposition (14.11) . 
Just as in (14. ip . to approximate || [X, rf(/9)]||2 we need to handle 

(27r) J 

. We will show that outside of some small ball around the origin, all phase space points are 

experiencing a uniform super-polynomial decay. 

Case 1: By (14.21) . for any e there exist a 6 such that |1| < 6 

-1(1 + e)(I|5-'-l) < ^'.(l) < -1(1 - 6)(l|i?-'-I), 

and 

supV^^(I) < sup-l(l-e)(l|S^'^l). 

|1|><5 |1|=<5 2 

Define the constant d, d = supm^a —(1 — e){\\B\). Define 5*5 j to be the set of phase space 

points (q, p) such that |q — sp| > | for at least a fraction of of intermediate times s in the 

interval [0,t]. Up to time t these points have a maximum decay factor of e'^^. It follows that 
for large times t these points have a super-polynomial and thus negligible contribution. On the 
other hand, points in S"* satisfy |q — sp| < 6 for all intermediary times s G [0,t] as long as 

3'* 

t > 4. This follows since the moving point q — sp requires a time interval of at least length 
t — to travel through an arc of B5{0). Hence in an additional time period of length ^/t, 

it can not travel the minimum distance y required to escape the 6 ball as long as t > 4. It 
follows that for all points in 5*1 and for all intermediary times s we have 

3'* 

-1(1 + e)(q - sp|S^'^(q - sp)) < " ^P) < "^(^ " '){^ " sp\B^^%ci - sp)). 

The region of points in Ss ^ is negligible so we have the asymptotic upper and lower bounds =F 
for our original expression as 

dprfqlq- tppe~-^o«'^(q-«p|^"'"(q-'^p))-/o'^^(iT€)(q-^p|«"'"(q-'*p)). 

By applying our results from (14. ip with A^'^ replaced by A^'^ -|- (1 ± e)B^'^ and letting e go to 
zero we get our asymptotics. 

Now we deal with the case where A^'^ is positive definite, but u is not assumed to have a 
density. Given any 6 the contribution from points in 5*4 will have a negligible effect on the 

2'* 

decay rate by the same argument as above through the term 

- / (is(q - sp|A'''^(q - sp)). 
Jo 
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Since points in Ss^ have that \q — sp\ < 5 for all intermediate times and by (I4.2p . for any e 
there is a 5 such that 2 ds'ip^{q — sp) is bounded above and below by 

^0 

By taking e less than the smallest eigenvalue of A^'^, we can apply (1) of (14. ip . and take the 
hmit as e goes to zero to get the asymptotics. 
The other cases are handled similarly. 

□ 

Now we give a theorem characterizing the asymptotic form of Tt{p). It doesn't require any 
new ideas as opposed to those appearing in the above proves. 

Theorem 4.4. Let the dynamics Tt satisfy the conditions of lji4-3{ ) for cases 1 or 3 and pt 
satisfy then 

\\UP)-Pth^^_ (4.6) 



\\Pt\\2 

Proof. By using flA.Sp . we can rewrite (14.60 using the quantum characteristic functions (pp^ and 
'pTt{p)- We can then apply the same techniques as in (14. 3 p to replace the Poisson noise by 
bounds using Gaussian noise. By Laplace methods, we can see that the numerator tends to 
zero faster than the denominator. □ 

Theorem 4.5. Let the dynamics Tt satisfy the conditions of \4-^ f^'^ cases 1 or 3 and pt{^,"v) 
be the probability distribution for the classical Levy process at time t, then 

(/dxdv|Wr,(p)(x,v)-p,(x,v)p)^ 

1 ^ u, 

(/ (ix(iv|pt(x, v)|2)2 

where >Vrt(p)(x, v) is the Wigner distribution for the state Tt{p). 

Proof. Taking the Fourier transform of the integrand in the numerator we are left with the 
difference between the quantum characteristic function of Tf (p) and the characteristic function 
of pt(x, v). The multiplication factor governing the noise in the formula for (y9r^(p)(q, p) causes 
points away from the origin to vanish. However, v^p(q + tp, p) is continuous and takes the value 
one at the origin, and this can be effectively replaced by the function that is 1 everywhere. □ 
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APPENDIX 



A The Quantum Characteristic Function 

The quantum characteristic function is defined as ¥'p(q, p) = Tr[lVq^p/9] for p G Ti(M'^). Weyl 
operators satisfy the multiphcation formula 

Formally, this formula follows from the Baker-Campbell-Hausdorf (BCH) formula. Using flA.ip 
with the characteristic function formula 

e5q-PTr[e^^-^e*P-^p]. (A.2) 

Since e***'^ acts as a translation operator by q in the x-basis, intuitively we can apply the 
formula for a trace to reach the equality 

e5q-PTr[e*''-^e*P-^p] = e-^q-P j dxe^P'Xx - q, x) . 
Now taking the Fourier transform in the p variable we get 

The Fourier transform of the q variable is by definition the Wigner distribution function 
Wp(x, v). Hence the quantum characteristic function and Wigner distribution function are 
related by a Fourier transform in both variables: 



^ j dq^ipe-^P-''-^^>,(q,p) = W,(x,v). 



(27r) 

Lemma A.l. Suppose p he a density operator, and J p G T2(-L^(M'^)) for 

G G {Xi, ■ ■ ■ , Xrf, Ki, - ■ ■ Kd\. Let qo, po G , with |qop + |poP = 1. Then we have that 

/i~^(W^h{qo,po) - I)P ^(qo ■ ^ + Po ■ K)p, and 

p(W^h{qo,po) - «P(qo ■ X + po- K), 

where the convergence is in the trace norm. 

Proof. Define the self-adjoint operator = qo ■ X + po ■ so we can write Wft,(q(,_pQ) = e^^^ . 
By our conditions on p and the triangle inequality, Hp is trace-class. Technically, Hp is defined 
as the bounded operator (traceclass even) determining the bilinear form B{g,f) = {Hg\pf), 
for g G D(if) and / G L^(M'^). In particular, the boundedness of the B implies that p maps 
arbitrary elements in L'^iW^) to D{H). 

Note that |/i~^(Vl^(/iqo,hpo)~-^)| ^ \H\- However, for two operators A, B such that < A < S, 
then pA^p < pB^p and < ||-Bp||i. To see that < let gj an orthonormal 

basis of eigenvectors for pB^p, then we have 

\\Aph < J2(^93\PA'P9,))'^ < J2((9j\pB'pg,)t^ = J2{9MpB'p)hj) = \\Bph, 
j j j 
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where the first inequahty above follows by writing pA^p in terms of its spectral decomposition 
and applying Jensen's inequality. Applying this fact with A = /i~H^(/iqo,/ipo) ~ ^) B = H, 
we have that 

\\h-\^'^''-I)ph<\\Hp\\,. 

Hence, {e^^^ — I) p is trace class. By the singular value decomposition, there exists a sequence 
of finite dimensional projections Pn such that HpPn converges to Hp in the trace norm. 

\\[h-\e'^'' -I)- iH)p\\, < \\{h-\e'^'' -I)- iH)pPr,h 

+ ||(/,-i(e^'^^ _ /) _ ^H)p{I - P„)lli (A.3) 

The second term is bounded by 2\\Hp{I — -Pn,)||i and we can pick a n large enough so that 
this term is smaller than |. On the other hand, the image of pP„ is finite dimensional and 
contained in the domain of H. Using Stone's Theorem [21] over that finite dimensional space, 
we can pick an h such that 

||(/,-i(e*^^^ _ /) _ iH)pP^\\^ < ^, and hence \\{h~\e'^'' -I)- iH)pPn\\oo < |. 
Hence we have the trace norm convergence 

^"^(W^h(qo,po) - iHp. 
Similarly p/i"HW^(/iqo,ftpo) - ^) ^ ipH- 

□ 

Lemma A. 2. Suppose p be as density operator and J p G Ti(L^(M'^)) for 

G G {-^1, ■ ■ ■ , Xd-, Ki, - ■ ■ Kd}. It follows that the first derivatives of <fp{ci, p) are bounded and 

continuous. Moreover, for qo, Po ^ we have the formula: 



(po) ■ ^('l'P)'^''('l'P) = ^'^{qo-/?+PO-X,p}(P'q)- 



Proof. Let |qoP + |poP = 1? h > 0, and W^(x,k) = e^^^~^^^^ be the Weyl operator for a translation 
by (x, k) in phase space. By the cyclicity of trace and action of Weyl operators 

V?p(q+ /iqo,p + hpo) = TT[W^q^p)W(^hcio,hpo)pW{hcio,hpo)]- 

We can write 

i((^p(q + hqo, p + /ipo) - </?p(q, p)) = TT[W(^ci,p)h'\W(hcio,hpo) - I)pW(hcio,hpo)] 

+ nW^^MWik^oAPo) - I)h-\ (A.4) 

By Lemma ^K^, /^"^(^^(feqo.ftpo) - I)p and piW^h^^^h^^^ - I)h-^ converge to i(qo ■ i? + po ■ X)p 
and ip(qo ■ K + po ■ X), respectively, in the 1-norm. Since VF(q,p) and W(^hcio,hpo) ai'e unitary, 
they are bounded in the operator norm, and the above expression converges to 

zTr[W^(q,p)(qo ■ i? + Po ■ X)p] + ^Tr[W^(q,p)p(qo ■ K + po • X)]. 

Since (qo ■ K + po- X)p and p(qo ■ K + po- X) are trace class, the expression above is bounded 
and continuous. Moreover, it can be written as 



qo 

Po 



^ ■ V(q,p)(/?p(q,p) = i^{^,.K+p,.x,p}iP,q.)- 



□ 
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Proposition A. 3. Consider the complex Hilbert Space H — L^(R'^). The map C sending trace- 
class operators p to their quantum characteristic functions (pp extends to an isometry from the 
Hilbert-Schmidt class operators T2(L2(M<^)) to L2(M<^ x M<^, j^dqdp). 

The above proposition can be proven be considering that L'^{W^) = r(C'^), where r(C'^) is 
the Fock space generated by C''. (pp is exphcitly computable in the case that p — \e{v)){e{u)\ 
where e{u) and e{v) are exponential vectors. It can then be shown that 

Tr[PiP2] = -^2^ J dx.p(pp^{x,p)(pp^{x,p), 

where pi and p2 are two non-orthogonal projections constructed with exponential vectors. This 
property can then be extended Hilbert-Schmidt operators in general. 
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